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Abstract 

We compute the Bethe equations of generalized Hubbard models, and study their 
thermodynamical limit. We argue how they can be connected to the ones found in the 
context of AdS/CFT correspondence, in particular with the so-called dressing phase 
problem. 

We also show how the models can be interpreted, in condensed matter physics, as 
integrable multi-leg Hubbard models. 



1 fomin@lapp.in2p3.fr, frappat@lapp.in2p3.fr, ragoucy@lapp.in2p3.fr 



LAPTH-010/10 
February 2010 



1 Introduction 



The nowadays widely studied Hubbard model was introduced in the sixties [HE] in 
relation with strongly correlated electron systems (see [3lH] and references therein for a review 
on the Hubbard model). In two and three dimensions, the Hubbard models are unfortunately 
not solved yet and few results have been obtained. In contrast, the one-dimensional Hubbard 
model was found to be integrable and its Hamiltonian was first diagonalized by means of the 
coordinate Bethe Ansatz by Lieb and Wu in 1968 [5]. Since then, there have been numerous 
studies on this model. 

Originally described as a model of spin up and down electrons on a one-dimensional 
lattice with hopping terms and nearest neighbors interactions (resp. kinetic and potential 
terms of the corresponding Hamiltonian), the Hubbard model has been then generalized in 
different ways. A particularly interesting method is based on Shastry's construction which 
was used to reveal the integrable structure of the Hubbard model [HIE]- The main idea is 
to couple two XX model i?-matrices through an interaction term depending on the coupling 
constant of the Hubbard potential. The proof of the Yang-Baxter equation for the obtained 
Hubbard R- matrix was given by Shiroishi and Wadati [8]. 

In this framework, a first generalization to the gl(n) case was proposed by Maassarani 
[Hj. The extension to the super algebraic case g/(n|m), mainly motivated by the appearance 
of the Hubbard model in the context of N = 4 Yang-Mills theories (see e.g. [ID]), was 
given in [11]. A general approach for deriving (super) Hubbard models was developed in 
|12j . The construction is based on the decomposition of an arbitrary vector space (possibly 
infinite dimensional) into a direct sum of two subspaces. The two corresponding orthogonal 
projectors allow one to define a i?-matrix of a universal XX model, and then of a Hubbard 
model using Shastry's trick. The QISM approach ensures the integrability of the models, 
and the properties of the obtained i?-matrices lead to local Hubbard-like Hamiltonians. 

We continue here the investigations of the one-dimensional integrable generalizations of 
the Hubbard model started in [12], where Hubbard-like models based on <?/(n|m) ©p/(n'|m') 
were introduced and the corresponding Hamiltonians explicitly constructed. The Bethe 
Ansatz equations were fully derived in the XX case, while only subsectors of the theory were 
examined in the Hubbard case, due to the complexity of the calculations. A step forward was 
accomplished in [13], where the case g/(n|m) ® gl(2) was investigated and the corresponding 
BAE determined. In this paper, we give the complete set of Bethe Ansatz equations for 
general Hubbard-like models. The structure of these BAEs is similar to the one of the usual 
Hubbard model, however with some (new) phases that depend on some Bethe roots that 
are quantized. These phases have to be compared to the ones introduced in the context of 
AdS/CFT correspondence [TU [15]. Hence, we believe that this construction may be seen 
as a first step in the construction of the integrable model underlying the super- Yang-Mills 
theories. 

Another interesting application of these models relies on their possible interpretation as 
multi-leg Hubbard models in condensed matter physics. This may be also a new way to 
tackle the problem of two-dimensional Hubbard models. 

The plan of the paper is as follows. In section [2J we remind the construction for general 
Hubbard models and we set the basic notations used in the paper. Then we present in section 
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[3] the set of Bethe Ansatz equations (BAEs) associated to these models, together with the 
energies and momenta. In order to make the presentation clearer, the organization and the 
details of the calculations, based on the coordinate Bethe Ansatz, are postponed in appendix 
[A1 The section S] deals with the thermodynamical limit. The Bethe equations in this limit 
are derived and the ground state energy is determined. Some applications of these models 
are presented in section [51 It is first emphasized that the corresponding Hamiltonians can be, 
after a Jordan-Wigner transformation, interpreted as multi-leg Hubbard models. Another 
point concerns the possible interpretation of the phases occuring in the considered Bethe 
equations in relation with the AdS/CFT correspondence. Finally, we conclude in section O 
on open problems. 



2 Description of the models 

Let us here remind the construction for general Hubbard-like models. For a given pair 
of algebras (gl(n), gl(vn)) there are roughly 



(n-l)(m-l)+l 
2 



non-trivial inequivalent models, 

corresponding to the choice of a projector in each of the two algebras. These choices are 
labelled by two integers (p, q), with 1 < p < n and 1 < q < m, that correspond to the ranks 
of the projectors. 

The Hamiltonian which we deal with in the calculations below is derived from a transfer 
matrix obtained from the standard procedure used for integrable systems. 

At first, we define the i?-matrix of an XX model based on the algebra gl(n): 

R 12 (\) = Ei 2 P 12 + S12 sin A + (I <g> I - E 12 ) P 12 cos A £ End(C n ) <g> End(C n ) (2.1) 
P 12 = J2 E i E 2 ( 2 - 2 ) 



P n 



'12 



^ ^ [EVEl' + El'E?), (2.3) 
i=i j=p+i 



where A £ C is the spectral parameter and p £ [1, n] is a free integer parameter that defines 
the model. E l J denote the elementary matrices (with entry 1 at row i and column j and 
elsewhere) acting in the x th copy of End(C n ). The XX- model i?-matrix (12.11) obeys the 
Yang-Baxter equation, is unitary and regular. 

The definition of the (generalized) Hubbard i?-matrix uses as basic ingredient the R- 
matrix of the XX model (or its generalization), which are coupled a la Shastry, the coupling 
constant being related to the potential u of the Hubbard model. The two underlying XX 
models may be based on two different algebras g/(n)t and gl(xn)i and depend on two different 
integers p and q [12]. Let us introduce the corresponding sets of integers 



A/" t = {l,2,...,p},A7 t = {p + l,...,n} and Af x = {1, 2, q} , = {q + 1, m} . (2.4) 
The i?-matrix of the Hubbard model based on the pair (<?/(n) , p ; gl(xn) , q) is given by 

it^Ai, A 2 ) = i4(A 12 ) 4a(Ai 2 ) + S -^r^\ tanh(/4) i&(Y 12 ) C n R l 12 (X' 12 ) C n (2.5) 
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where A12 = Ai — A 2 , A' 12 = Ai + A 2 and we introduced the diagonal matrix C a (a =7 or |): 

c* = E E i j - E P « ( 2 - 6 ) 

The coupling /i' 12 = /i(Ai) + /i(A 2 ) is based on the function h(X) such that 

sinh(2/i) = usin(2A) . (2.7) 

The i?-matrix (I2.5P is symmetric, regular and satisfies the unitary relation. Moreover, 
when the relation (12. 7p holds, the i?-matrix (12. 5 p satisfies the Yang-Baxter equation: 

-Ri 2 (Ai, A 2 ) i?i3(Ai, A 3 ) -R 2 3(A 2 , A 3 ) = Rli{\-2, A 3 ) R\i(Xi, A 3 ) -Ri 2 (Ai, A 2 ) . (2.8) 

Being equipped with an i?-matrix with all required properties, we can proceed to define the 
corresponding quantum integrable system, by performing the following steps: monodromy 
matrix, transfer matrix and Hamiltonian. The L-site monodromy matrix is given 

T a<bl ... bL> (\) = Rli(X, 0) . . . Rli L (X, 0) (2.9) 

and its transfer matrix is the (super)trace in the auxiliary space: 

t(X) = tr a T a<bl ... bL> {X) . (2.10) 

Then the generalized Hubbard Hamiltonian reads 

d 



— E Hx,x+1 (2.11" 
A=0 x=1 



with ^ 

H x ,x+i — (^P)t x,x+i + (^P)ix,x+i + t C^ x Ci x , (2-12) 

where we have used periodic boundary conditions. 

We consider generalized (gl(n),p ; gZ(m),q) models, containing thus four different types 
of particles n 7, tx 7, n I and n I, each type being 'colored': the n 7-particles have 'colors' 
1 t, 2 |, P T; the 7f t-particles have 'colors' (p + 1) 7, n t, while the 'colors' for ir I and 
7f ^-particles are 1 J,, q I and (q + 1) |, ...,m j respectively. 

The L-site Hamiltonian is given by 

L 

H = H (( SP )t + ( SP )^ -.-+1 + \ C ^ . mod L ( 2 - 13 ) 

x=l 

where u > and 

( N -/')- ,,., = E E { E IHU + > ( 2 - 14 ) 

= e E { E t E t^ + C^C+i) > (2-i5) 



The corresponding system will be called a (gl(n), p; gl(m), q) model, or a ; ™J -model. 

Note that the ^ n " p ; J -model is equivalent to the ; ™ J -model. The usual Hubbard 
model is the ^)-model. The models introduced by Maassarani j9] are the (" ; ™)-ones. 
Generalizations to the superalgebra CclSG clS done in [12] will be noted ; '^j') -models 

with obvious notation. The generalization given in [II] corresponds to (™f ; m ' 1 ;') -models. 



The symmetry algebra of the ( £ ; m )-model is a ©g/(n — p) ffig/(q) ©g/(m — q) algebra. 



Remark that one can add to the Hamiltonian chemical potentials 

L 

E fa E ^ + ^ E E t + k E . + At E ^ 

x=i ieA^ j- e ^ ie/V t j 6A f t 

without perturbing integrability nor symmetry. 



3 Bethe equations of the ; q J -model 

We present in this section the main result of the paper, namely the Bethe Ansatz Equa- 
tions of the model under consideration, see Hamiltonian f)2.13p . The details of the compu- 
tation are postponed in the appendix and can be skipped in a first reading of the paper. 

The spectrum of the generalized ; ™j -Hubbard model is 

E = ^(L-2N*) + 2 Yl cos ^) C 3 - 1 ) 

where the Bethe roots ki are parameters solution of the Bethe Ansatz equations (see below) 
and L is the number of sites. To present them, we introduce integers < K < N < L, 
Nfa, iV-j-s., N± n , > such that K = N^ n + + N± n + N^, and sets of integers 

= {1, • • • , N t7T } , = {N tn + 1,...,N U + N n } , (3.2) 

M n = {N^ + N n + 1,...,N U + N n + N u }, (3.3) 

M n = {N u + N n + N l7T + l,...,K}; (3.4) 
A = U A w i = {ax, a 2 , • • • , a^} U {aN u +N n +i, o-N u +N n +2, • • • , a N^+N n +N^} (3.5) 

The integers a>j are such that Oj ^ cij for i G M^, j GM4 and 

1 < ai < a 2 < ... < a Nu < N , 1 < a^ + jv t# +i < ... < a Nu+Nn+Nin < N (3.6) 

Then, Bethe Ansatz equations are 
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9 ,- m 1~7r 

e ifc,(L-iv ts ) = ( _ 1 )Ar t .-i e ^ mt7r = i ; ... ;A r^ for jGA?rt! (3.7) 



e 



- (-i)^^- 1 n e ~ tkam n isin ! j t^ + ! n ^ > ( 3 - 9 ) 

- LJ - 11 zsinfc,- +iAi - 7 - LJ - 



z sin A;,- + ia 



for j€[l,iV]\ 

2tN vl JV 4-* -J, 4.(m-l) 



h N^N lm = eN -T.^ % iGiM ^ a ^(6 z )<ar^(M (3.10) 

with 1 < m < ... < rtN^-N^^ < N)*, 

II a ,. A| - = n 83in<: , + iAm ^ . for (3.11) 



with 1 < ni < ... < n N ^- Nu < N n 

We have chosen e^^e 1 ^ as the reference state at the first level of the Bethe ansatz (see ap- 
pendix . so that all states with e 2 ^,...,e"^ and e 2 ^,...,e m ^ appear as excitations ("particles") 
above the reference state. N is total number of such particles. 

The parameter N-fi is the number of % j" particles for % = 2, n and is the number 
of j I particles for j = 3, m. = Y^=2 an< ^ = ^21=2 coun t the number of 
7r-particles with spin up and spin down respectively. Nfa = Y^=p+i ^t* coun ts the number 
of spin up 7f-particles, while = Y17=q+i counts the number of spin down 7f-particles 
that are not of type m J,. The reason for this latter choice will become clear in the following. 
In the same way, K = N^ n + + N± n + is the total number of particles that are not 
of type m |. 

For given integers N^ n , N^, and N^, the phases bj and the integers rij and nj 
correspond to the different "colors" that can have particles of a given type (t tt, t ^> 4 n or 
i 7f types). 

The integers a,j (entering in the sets A^-j- and A^) define the order between momenta of 
the t Ti" and J, tt particles. This order is preserved (up to a cyclic permutation) by the action 
of the Hamiltonian on the wavefunction. 
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4 Thermo dynamical limit 



The Bethe equations given in section [3] differ from the (usual) Hubbard model's ones by 
some phases. In this section, we study them in more detail and look for their thermodynam- 
ical limit L — > oo. 



4.1 Simplification of Bethe equations 

The Bethe equations f l3.7j) -( T5TTTj) can be partly solved for Bethe roots kj with j e A^-j-UA^ 
and bi with / G Ai^i'- 



2vr / 




l - jv t A 




m tn = 1, .. 








i - ^ - 




m u = 1, .. 





'I- < ... < 4; 



F-4.7T T^T 
"I < < ^ 



1 < ni < ... < < JV W I 1 < If < ... < /£, < Af w + iV ta 

^infr-HA,,,- = C II iAm -iA,- ' fOT m £ ^ < 45 > 



where the phases $ and \& are defined by 



2 N Xn J L-N^-N^ 

-y — ^ + v — ^ — + y — - — (4-6) 



and 



^ + N i* + _ % fj ^_ (A7) 



with 1 < m < ... < n Nf ^ Ntn < N> 
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We recall that 

A = U A^ = {a 1; a 2 , . . . , a<v t7r } U {a/v^+jv^+i, a iv t7r +iv tif +2, • • • , a N tn +N n +N in } (4-8) 
where the integers are ordered according to the inequalities 

1 < ai < a 2 < ••• < ciAr t7r < N and 1 < a Nu+N ^ + i < a Nu+N ^ +2 < ... < a N ^ +Nn+Nu < N 

(4.9) 

and Oj 7^ for i G M.^, j G .M^. 

In order to coincide with the results in [16], we now slightly change the notations 



Aj -> -Aj and u -> -4C7 (4.10) 
as well as the sign of the Hamiltonian for the energy to be equal to 

E = —U(L — 2N^) — 2 cos (^) ( 4 - n ) 

leM^uM^i 

The ground state of the model is given by the real values of the Bethe roots kj and Aj. 
Following Lieb and Wu [5], we take the logarithm of the Bethe equations (I4.4p and (J475]): 



k 3 L = 2^ + J i ) + ^g( A< J infcj ), je[l,N)\A (4.12) 



N A . , N n 



Aj — sin kj 
U 

i=l 

±e(^p^-) = 2, W + *) + f>(^-^), ieMn (4.13) 

3=1 k=l 

where #(2) = 2 arctan(x) el — 7r, 7r] and we take the cut off for the logarithm as - log ( 

i \x — i, 

n-6(x) G [0,2tt[. 

The quantum number Ij is integer or half-integer depending whether + + N± m — 1 
is even or odd, similarly Jj is integer (half- integer) for + 1 even (odd). 

4.2 Thermodynamical limit 

N - N 



We consider the thermodynamic limit L — > oc where the particle densities 



f £T7T - 1 CT7T 



are kept fixed (a =t, J,). Considering the phases $ and \& in the thermodynamical limit the 
$ $ 

ratios — and — do not vanish and depend on the particle densities. 
L L 

In this limit, the real numbers kj and Aj are close one to each other: kj+i — kj — > 0, 
Aj + i-Aj -»■ with L ->■ oo. They are distributed between -Q($) and Q( < &) < 7r and —B(ty) 
and _B(\I/) < oo for some Q(<&) and B(ty). In the small intervals dk and dA, the numbers of 
kj and Aj are Lp(k)dk and La(A)dA respectively, where p(k) and cr(A) are density functions 
to be determined. They are normalized as follows: 

p(k)dk = - 7 ^ and / <r(A)dA = ^ (4.14) 
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The counting functions I{k) and J (A) are defined as usual from Ij and J, in the continuum 
limit: 

1 f Q( -^ /A — sin h\ 1 r B w /A — A'\ 

J(A) = f / ^)<^F) - * - f / rfAV(AO^^-A) (4.16) 

2tt y_Q ( $) \ U J 2ix J_b W \ 2U J 

These functions are also such that I(kj) = Ij and J(Aj) = Jj. 

Since J(A; + dk) — I(k) counts the number of k values between k and k + dk, we have 
= Lp(k), and similarly, — - — — = La (A). 

Now taking the derivatives of (14.131) . and considering the fact that the phases $ and ^> do 
not depend on the Bethe roots k or A, we get the equations on densities, which are the same 
as Lieb and Wu in [5]: 



cos A; f , . . „. /A — sin k\ 
2np(k) jj J dAa(A)6'( J (4.17) 



where k E [-Q,Q], A e [-B,B] and 9'{x) = d9/dx(x). 
4.3 Ground state 

We consider the "half-filled band" limit, defined as 

N-N u -N lw = L and 2N n = N - N u - N l7T . (4.19) 

Using the same arguments as in [5], we obtain Q(<&) = vr and B(ty) = oo. This can be seen 
as follows. Taking the normalization relations for p(k) and a (A) (I4.14p and inserting there 
the relations between the counting functions and the densities, we get: 

I(Q)-I(-Q) = N-N u -N ln , J(B) - J(-B) = N n (4.20) 

Inserting the definitions of the counting functions (14.161) in these equations, the following 
conditions arise: if Q — > 7T then N—N^—N^ = L and if B — > oo then 2N^ = N—N^—N^. 
This limit allows us to find the solution for the densities p(k) and a (A) by Fourier transform: 

a (A) = -L f X 'dw j o(w ) A e [-00,00] ( 4 .21) 

Z7T Jo COS{WU ) 

1 cos/c f°° , cosfw sin /c) T . . , r , . 
Po(k) = - + — j( rf W g2 L + ! J °M , fe e [-7T, 7T] (4.22) 

i r 00 

with zeroth order Bessel function Jq(x) = — dw cos(xsmw). 

7T Jo 

The ground state energy is then equal: 

£ = -U{L-2N*)-2 V cos(fc,) = -U{L-2N«)-AL f°° dw ^ffi^ ( 4 - 23 ) 

L — ' / n w\e AVW + 1) 



X f°° 

with order one Bessel function J\{x) = — dw cos(xsinw) cos 2 w. 

7T Jo 

4.4 String hypothesis 

From the study of the ground state, it is tempting to conjecture that in the thermody- 
namical limit, the Bethe parameters line up into so-called strings, as for the usual Hubbard 
model. We remind that the string hypothesis states that all regular solutions {kj} and Aj 
of Bethe equations (13.7jl - fl3.lip consist of three kinds of configurations: 

1. real kj <G [— n, n]; 

2. A's combined in A-strings (of arbitrary length n): 

= A n a + iU(n + l-2j), j = l,...,n, (4.24) 

3. 2n /c's and n A's combined in k — A-strings: 

A£ = A' a n + iU(n + l-2j), j = l,...,n, (4.25) 

(it -&rcsm(A' a n + iU(n + l - j)), if j = 1, 3, 2n - 1 
k n J = \ &rcsm(A' a n + iU(n - j)), if j = 2, 4, ...2n - 2 



i \ a \ j . 

I 7r — arcsin(A' a n — inU) 



if j = 2n 



The string centers A", A'J 1 are real, the parameter n indicates the length of the string and 
a priori goes from 1 to infinity, the parameter a counts the strings with the same length: 
a = 1, M n , M n being the number of strings of length n. The branch of arcsin(x) is chosen 
between —ir/2 and n/2. 

Using the string hypothesis inside the Bethe equations, especially for equations (14.41) 
and (14. 5p . and taking the logarithms, we arrive at the following form of the Bethe Ansatz 
equations for the (real) centers of the strings, called discrete Takahashi equations: 

M , , M' 



n=l \ a=l a=l 

N—2M' . In _ . , 

2L^(arcsin(A' n + inU)) = 2TrJ' a n - 2vr($ + V) + ^ fl( a w ^ 



no M' 

°° m , A in \'rn x 

+ (427) 

m=l 6=1 



N ~ 2M ' /A» aini-.N °° Mm / A n A m \ 



i =1 m=l 6=1 
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where j e [1, N - 2M'\ \ A in fl4T26j) and M' = Yln=i nM n- The cut-off for the logarithm is 
taken as in section I4TT1 The function Q nm (x) is defined as: 



m-\-n -I 

2 



<W*>= £ «|)^ + "(5(iT3))) < 4 - 29 > 

where <5^o is when / is 0, and 1 otherwise. 

Looking at the range of the functions entering the BAEs, we can give bounds for the 
integers I jt J™, J' a n : 



L 



-.CO j . 'SO 

- J2( M m + M'J - min($) < /, < - - - £(M m + M'J - max($) (4.30) 



2 2 ^ y m/ y ' ~ J ~ 2 2 

m=l m=l 
T/n ^ T/n ^ T/n 



JZn<J'a n <Ca X with (4.31) 



min 



■1 



max 



-UL-N + 2M' + N u + N u -l-jr M' m t nm j + max($) + max(^) 

\ m=l / 

^ ( L - A + 2M' + iV t7r + A^ - 1 - ^ M^t nm ) + min($) + min(W) 

\ m=l J 



2 

J n ^ n <J n a <J n max with (4.32) 



jn 

mm 



J, 



max 



~(n- 2M' - N u -N u -1- M' m t n ^ + max(^) 
\{n~ 2M' - A t?r - ^ - 1 - MLt nn )j + min(tt) 



with t nm = 2min(n, m) — 5 nm , min(v4) (resp. max(v4)) is the minimal value of A and 
Ij, J™, J'™ are such that: Ij is integer (half-integer) if 5^^ > =1 (M m + M' m ) is even (odd); J^ n 
is integer (half-integer) if L — (A — 2M' — N^ n — N± n ) + M' n is odd (even) and J" is integer 
(half-integer) if (A - 2M' - N u - N u ) - £~ =1 M ^mn is odd (even). 
The minimum and maximum of $ and \P are: 

max( $ ) = ^±3^1 N l + l +N N»-N«~-i)(N«~i) + W 

mm <P = — A 1^ 

1 ' 2 L — A^j — Aj, m in 

+ 2(A„ + A, m ) (4 - 33) 
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N n + N ijf + N im , fAT AT AT , fAT AT .. , N u 



maxW = ■ * ^ ^ + (7v t7r + iV ;7r + iV^ + (7V tff -iV tw )) + 



N u (N i7r - 1) + 1 (iV> - N im ^)(N^ - N im ^ + 1) - N^(N^ - 1) 



(Ar t --iV tn )(iV t# -iV tn -l) 



2N n 

{ 1 2 L-iV t# 2 (N ijt + N lm ) 



(4.34) 



| (iV t ,-Ar tn )(Ar t ,-iV tn + i) 

If the string hypothesis is still valid, these equations should be the ones to consider for the 
study of excited states above the ground state described in section 14.31 The first numerical 
studies we have performed seems to indicate that when the phases \l/ and <3> are not zero, 
the string hypothesis needs to be modified: the spaces between the imaginary parts of the 
Bethe roots seem not to be integers anymore. However, more detailed numerical studies are 
needed to confirm it. If this alteration is confirmed, than the above Bethe equations will be 
valid only in the " zero phase subsector" . 

5 Applications 

5.1 Multi-leg Hubbard models 

In this section we consider a particular ; j^J Hubbard model. The Hamiltonian 
(I2.13P for superalgebras as introduced in [12] can be written as 

L „ 

H = -J2 ((£P) t „*.+i + (EPW+i - J C^C^J, mod L (5.1) 

2 = 1 

with the gradation introduced in kinetic terms: 

(XP) ax , x+1 = + with a=t,| 

J=2,3 

We use the Jordan- Wigner transformation (see [IT] or [IB] for more details): 



El\ = (1 - < t Xt Z H> E Z = C 1 - <Hi ( 5 - 2 ) 

El\ = -(1 - nf t )4 t Z^, E 3 a \ = -(1 - n^)4 4 (5.3) 

with < CT = 4 CT c aCT , n d au = dlJ aa and = Utii 1 ~ 2 ™u)- 

The Hamiltonian ( 15. ip finally can be presented in the following form: 
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H = H Hub {c,c\U) + H Hub {d,d\U) + V mt {c,c\d,rf,U) (5.4) 

L L 

H Hub (c,c\U) = -^^(4 Cfe+1(7 + c t +iCTCfeCT ) + ^^(i -2^(1-2^) (5.5) 

i=l cr=t,l i=l 



with 



V int {c,c\d,d\U) = ~E E ((4a c k+ia + cl +la c ka )(nl a n d k+la - n d ka - 



n k+l,a) 



i=l a=t,i 



(d ka d k+la + d\ +la d kcT )(n c kta n c k+1 ^ 



n k.p n k+l, 



u L 

- E ((1 - 2nJ t )(l - 2<) + (1 - 2<)(1 - 2^ 



■ . 

fc=l 

2(1 - n* t - <)(1 + 2n^<) - 2(1 + 2n c fcf < t )(l - n 



+ (l + 2< t nft)( 1 + 2 <;0) ( 5 - 6 ) 

The Hamiltonian ( 15.4ft can be interpreted as two periodic lines on which the electrons (de- 
scribed by c, for the first line, and by d, d> for the second one) interact via the usual 
Hubbard Hamitonian (15. 5p . plus a term of interaction Vi nt between the two lines. Thus, one 

gets a 'double-row' Hubbard Hamiltonian, and considering more involved ; ™J -models, 

one can construct multileg-Hubbards models. In this way, we construct an almost two- 
dimensional Hubbard model that is still integrable. 

The eigenfunctions for this Hamiltonian are made of creator operators cl,ct,dl and d\. 
They can be written in the following form and correspond to the solutions constructed in 
appendix [AJ 

*N M , Ni = £ E E £ *v (i u ( v (3 u (4) ) 

SW £(2) 2(3) iW 

k *~ I k ' I ' i i *~ j i j ' I 

Ni N 2 N 3 N 4 

ih^ n< ( 3) t ih^ i° >= 

i=l J i=l 3 j=l 3 3=1 3 

TVi N 2 N s N 4 

and A = (3tT^3t, uT^H, 2T~^2~t uT^H) ■ (5.8) 

where <f)[A] is the eigenfunction given in (1A.1I) . The particles c|,c|, dj and dj correspond to 
3 t, 3 4, 2 t and 2 j. respectively. The particles c£ and d£ can be identified with a (spin up 
and down) electrons. 

The energy of the excited state $Ni,N2,n 3 ,N4 reads 

N 

E = u(L- 2N) - 2 E cos h with N = N 1 + N 2 + N 3 + N 4: 

i=i 
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where the parameters k\ are Bethe roots defined by equations: 



e ik jL . i-r A/ — sin ki - lU 



nAl — bill K n — LU -I— r , „ r ^ _ Tl 

J—J-^Jlb,, for 3 e[l,N] (5.9) 



1=1 J 1=1 



6f 2+7V4 = 1, arg{bi) < arg(b i+1 ), with i e [1,N 4 ], (5.10) 

n\i — sm ki — lU v-jv 3 = . -i — |- -i — r A; — A m — zz(7 

j=i 1 j i=i ™=i 1 m 

with 1 < ni < ... < n Ns < Ni + N 3 , for I e [l,Ni + N 3 ] 



5.2 Comparison with AdS/CFT 

The (energy) spectrum of the (usual) Hubbard model has been shown [10] to reproduce 
correctly the spectrum of the dilation operator in the su(2) sub-sector up to three loops. 
The perturbation theory is done at g — > on the SYM-side, while it is done at u — > oo on 
the Hubbard's side. Starting at fourth loop, the corrections [2] differ with the Hubbard 
model. The string Bethe equations (for the su(2) subsector) read 

/„(„, i i\\ L m . | m = |cotan(ip,) 

x Jpi±4) .n ^-^' e^) with / / 2 r~^\ (5-1-) 

where 9(uk,Uj) is the so-called dressing phase that takes care of the discrepency to the 
Hubbard model. In general, this phase has two different origins. One is just related to the 
asymptotic Bethe ansatz and the scattering matrix of the model, and the other cause comes 
from the so-called wrapping problem, occuring when the order of the expansion in g is bigger 
than the length of the chain, L. When considering the su(2) sector and the Hubbard model, 
the wrapping problem has been already 'solved': the Hamiltonian is local, and it is the 
expansion in term of the coupling constant that makes appear different spin chain models 
with long range interactions. Hence, for this sector, it is only the dressing factor that has to 
be implemented. 

In the Bethe equation (13. lip , such a phase occurs, that depends on some Bethe param- 
eters. Unfortunately, for the present models, the Bethe parameters involved in the phase 
seem not to be of the type Uj. However, we believe that the present construction could be 
a first step for the construction of an integrable model possessing such a dependence. In 
particular, since the models we construct can be interpreted as multi-leg Hubbard models 
(see section |5TT|) . they are free from any wrapping problem. In the context of SYM theories, 
one should stick to the single circle interpretation, with electrons possessing some internal 
degrees of freedom (as it is done in the Bethe ansatz of appendix [X]). Then, to get a true 
su(2) model, one would have to integrate over the internal degrees of freedom, to get an 
effective model. 

To be more appealing, we reformulate the Bethe equations (I5.9p - (15.1ip with the notation 
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and half- filling constraint: 



2U = -^= ; ^ = u t ; N = N 1 + N 2 + N 3 + N A = L ; N x + N 3 = M (5.13) 
gy/2 2.U 

Forgetting for a while the state multiplicity, eq. (I5.10p . we also set 

Y[bi = e^ e2i and e 7 ^^ E ^ fi< = e^ 13 (5.14) 



t=i 



Then, we get: 



jHue,, rr«i- fl 1 ^ sin fc.- 



e ^ = e ^^TT ^ yv ; MUft » ? for j = l,...,L (5.15) 
_^ — gv2 sin kj + | 



— p\/2 sin 



TT - " r J - = e^"*^*" TT for Z = 1,..,M (5.16) 

M u, - SV^ sin ^ + f n ; - Um + z 



The first equation can be used to determine the the Bethe roots kj, while the second one 
has to be compared with ( 15.121) . 



6 Conclusion 

We have presented the Bethe equations for generalized Hubbard models, using the coor- 
dinate Bethe ansatz. They look similar to the usual Hubbard model ones, but a phase that 
appears in the equations. This phase has to be compared to the one appearing in the string 
Bethe equations. Unfortunately, the present phase is constant (more precisely depends only 
on the non dynamical parameters associated to the n l particles,), while the one appearing 
in the string context depends on the Bethe parameters (associated here to the 7r t particles). 
We believe our construction can be used as a first step for constructing an integrable model 
that would reproduce the expected phase. An open problem is thus to look for an amend- 
ment of the construction to provide /c-dependent phases for AdS/CFT. In that respect, let 
us note the studies done in [TT] that provides an alternative way to get generalized Hubbard 
models. 

The thermodynamical limit has been also discussed and a string hypothesis presented. 
In some cases (zero phase subsectors), it seems clear that the string hypothesis is correct 
and similar to the one of the (usual) Hubbard model. However, when the phases are not 
zero there are presently not enough data to conclude on the validity of the string hypothesis. 
Thus, we think that more numerical studies are needed to confirm or modify the string 
hypothesis in the thermodynamical limit. 

Finally, we think that our models can be used in condensed matter physics to define of 
multi-leg Hubbard models. Contrarily to the real two-dimensionnal Hubbard model, these 
multi-leg Hubbard models are still integrable. Hence, they can be used to get more insight 
on the two-dimensionnal model using integrable methods. 
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A Coordinate Bethe ansatz 



The derivation of the Bethe Ansatz equations of the model is based on the use of the 
coordinate Bethe Ansatz. The diagonalization of the Hamiltonian (12.1 3ft . which involves 
n + m types of particles, is done in two main steps. The Hamiltonian eigenfunctions are 
written as linear combinations of plane waves, whose coefficients are found to be solutions 
of a new eigenvalue problem, the Hamiltonian of which shows up as a chain of S'-matrices. 
This is the first auxiliary problem. 

In order to diagonalize this auxiliary Hamiltonian, the corresponding wave function is 
decomposed as excitations above a suitably chosen pseudo-vacuum with some coefficients. 
Using a recursive representation of the auxiliary Hamiltonian, recursive relations between 
these coefficients can be found. The resolution of these relations depend on the effective 
structure of the Bethe roots, since the choice of the pseudo-vacuum breaks the symmetry 
between spin up and spin down particles. The different cases are then considered. 

A.l Coordinate Bethe Ansatz, level one 

The states corresponding to N excitations are written as 

N 

x6[l,L] i=l 

with A = (2 |, 2 f , , n t, • ti f , 2 |, .?., 2 |, m ~, .?., m |). The sum is done on x = 
(xi,X2, ...,Xn) without points where two particles with the same spin coincide. The vacuum 
is chosen as 

0o=ri e ^- ( A - 2 ) 

k=l 

Acting with the Hamiltonian ( 12 . 1 3[) on the excited state, we get the eigenvalue equation 
for the \& function: 

(*[ x - e i + (^e^A- + [x + e, - 5+e n ,A]A+) + 
(J(L-2JV,)+u J Yl ^=^n)-^)^[x,Z]A 3 = (A.3) 

where e m is an elementary vector in with 1 on the m position and elsewhere. The 
notation \l/[x — e; + ^e,,.] means that, if under the Hamiltonian's action some 7f-particle 
(at position %\ with some spin a{) is moved to a place already occupied by some other 7r- 
particle with same spin, then they exchange their places. The symbols Af, A 3 represent the 
exclusion principle for 7f-particles: they are equal to zero if two such particles with the same 
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spin coincide on the same site. Explicitly, they read: 

a - = n n 5%{ - xi * x «) n st ^ ± ^ ^ ± i) , ca.4) 

l^m n^m I 

A 3 = Y[5H Xl ^x n ), (A.5) 



where 



8$(xi 7^ x n ) = 1 — 8{xi — x n ) for l,n E M.^ or M.^. (A. 6) 

We assume the Bethe hypothesis for the general solution of ^[x, A]. Dividing the coordinate 
space (xi, x 2 , .., xat) into AM sectors, we write for xq(i) < £q(2) < •• < ^Q(jv)> 

*g"[x,3|= Yl ^{PQiP- 1 )^ , P = P^P' (A.7) 

P — Pit Pit 

where k±, .., are unequal numbers (the Bethe roots), Q = [<?i, <?2, ••, <Ziv] is an element 
of the permutation group &n and Pkx = J2i^p(i) x i- We consider the permutation P in 
its factorized form P = P n ^P 7T P ji ., where P w ^ is a global and fixed permutation of tt and tt 
particles, while the terms P n and P^ permute tt and 7f particles separately. 
The energy depends on the "global" P^ permutation: 

E p " = 2 Yl cos(k PMl) ) + ^(L-2N^). (A.8) 

The coefficients <&j(PQ, P _1 ) in ( 1A.7I) are not all independent. Indeed, using the symmetry 
of the wave function and the application of the Hamiltonian represented in (IA.3j) . it is possible 
to reduce their number in several cases: 

1. For identical particles of any kind (tt or tt), the wave function satisfies the following 
symmetry property 

t x > 3) = *Sn ii+1 [nQ (i )Q(i +1) x, n Q(0Q(i+1) A) (A.9) 

that leads to the relation 

<^(PQ, P^ 1 ) = <^(PQ, (PDgwgfi+i))" 1 ) ( A -!0) 

where II a b is the permutation of objects a and b, the indices Q(i),Q(i + 1) correspond to 
identical particles. 

2. For particles with the same spin we also impose some kind of weak "exclusion principle" 
such that the wave function vanishes if two particles with the same spin coincide (but the 
particles can differ by their tt or tt type). This principle is verified for 7r-particles if we use 
( 1A.3I) with two particles at positions xq^ = £q(i+i) — 1 (in the sector £q(i) <C .. <C XQ(i) < 
XQ(i+\) <C xq(N), the notation <C means that the particles are far enough from each other). 
We generalize this condition to any kind of particles. It leads to: 

*q 1 x ^k (l ,-o ( , + i) = > A Q({): A Q(i+1) GA/juJV; or AT t u AT 4 (A.ll) 
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that is to say 

%(Pg, P- 1 ) = -^(PQU ll+1 , (PIL Q{i)Q(i+1) )- x ) (A.12) 

for any i such that Aqq, Aq^+i) G A/^ U A/4, or ATf U ATj,. 

3. For particles with different spins, there is no exclusion principle and they can be at 
the same site. Thus, let us consider the case: £q(i) <C •••• <C x Q(i) — x Q(i+i) ^q(at) and ask 
for the continuity of the wave function on the boundary of the two domains Q and QUu + i: 

^"[x, A] XQ(t)=XQ(i+1) = Wgn^Jx, A\ Xq( . )=Xq( . +1) . (A.13) 
For two 7r-particles, it implies the relation, for Aqu) G A7*j-m and Aqa +1 ) G A/"^), 

$ 3 (pg, p- 1 ) + %(pgn ii+1 , (piiq^+d)- 1 ) = 

%(pg, (png^Q^D)- 1 ) + %(pgrr m , p- 1 ) (A.i4) 

while for two 7r-particles it leads, for Aq^ G A/^) and Aq(j +1 ) G A/^(^) , to 

%(pgn ii+1 , p- 1 ) = $^(pg, p- 1 ) . (A.15) 

When one particle is of type tt and the other one of type tt, we find, for Aqu\ G A/^-m and 
A Q { l+ i) g AT^t), 

$ z (pgn ii+1 , p- 1 ) = $^(pg, p- 1 ) . (A.ie) 

There is another condition for 7r-particles when using (1A.3[) . This is the usual relation on the 
coefficient ^(Pg, P" 1 ) obtained in the (? ; ?) Hubbard model for spin up and spin down 
interacting electrons: 

(*q * [x - e m , A] + [x + e m , A]) + - 2JV ff ) + u - P) [x, A] + 

Skipping the intermediate calculations and combining the results with flA.14[) . we can write 
the conditions on $(Pg,P _1 ) in a matrix form: 



(A.17) 



^(rW^p- 1 ) \ = ft ab r ab \f ^{pq.p- 1 ) 
$(n a6 pg,(n afe p)- 1 ) y ' V r ^ t ab J \<s>(pQ,(n ab P)- r 

with a = PQ(i), b = PQ(i + 1) and 

2i(X a -X b ) -u 
u + 2t(X a -X b ) u + 2t(X a -X b ) 

These equations hold for any type of 7f-excitations (any value of Aq.,Aq. +1 being 7r-particles) . 

4. Finally we consider the interaction between w and tt particles with the same spin. 
Let Aq(i) be a 7f-particle, Aq^+x} be a 7r-particle, with coordinates xq^ = XQ(i+i) — 1, thus 
5f = 1 and 5f = 0. Using the equation (1A.3j) we derive the relation: 

*§L + J X + e W) " e Q(m),^] - *Q *( x + e Q(i)-^ = (A. 19) 
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that implies the following condition on ^^(PQ,P x ): 

^PQ(i)PQ( i+ i)PQ, P~ l ) = e***w*j{PQ, P^ 1 ) . (A.20) 

We rewrite now all obtained conditions fjAlU]) . (Rl2)l - (j07)l and fjAT20|) m a more 
compact form. Let us introduce, for P' = PQ G &n and Q' = P _1 G 6jv, 

A 7 " 

$(po = x>*^w) ik q ' w (A.21) 

»=1 



where the summation is over all types of excitations and all corresponding sectors. The vector 

[\f = i &i Q ' w belongs to V± <g> ... £g> Vn, where V is spanned by {e 2 ''', e 31 ", e n ^, e 2 ^, e 3 ^, e m ^} 
and represents one type of N excitations. The ordering of the particles is chosen such that, 
for Q' = id, the vector YliLi e t* * s taken as in flA.lj) : 



iV 



J ef 4 = e 2t ® ... <g> e 2t ® e nt ® ... <g> e nt <g> e n ® ... ® e 2; <g>... ® e m4 ® ... <g> e" 4 (A. 22) 

t=l 

To clarify the notation for $(P'), we give as an example the case N = 2: 

n m 

HP') = J2^rt)(P'M)eM + J2^dP'^)eM 

3=2 3=1 



+ E ( P '> ^ ef ef + (P', n 12 ) ef 

2<i<fc<n 

+ E ( $ W.*U ( P '> ef ef + $ WiJ4) (P', n 12 ) ef 

2<i<fc<m 
n m 

+ E E {^M)(P', ^) ef ef + * 0tM) (P / I H la ) ef ef ) (A.23) 

j=2 fc=2 

Then, all the relations can be expressed in a matrix form. For N = 2, it reads 

$(Il 12 P) = S { 1 l 2 \k 1 ,k 2 )$(P) (A.24) 

where S^(k\, k 2 ) (simply denoted by in the following) acts on elementary vectors e^ef 2 
as 

1) for 7r-particles: 

s$ 4'4 2 = -4 2 4\ for a u a 2 g Af m 

SS = for A x G A%), A 2 g A^ (t) (A.25) 



2) for 7r-particles: 

ii _Ai _Ai r T a s- T, 

Hi) 



S$ ef'ej 2 = -ef 2 ej\ for A X ,A 2 G AT t 



S{ 2 } ef ^ = t 12 ef ^f 2 + r 12 ef 2 ef for A, G A7 m) , A 2 G A7 4(t) (A.26) 
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3) for mixed it and 7r-particles: 

Sg> e^ep = e-^e^ep and Sg> epep = e ik *e?*e$\ for A, e A% } , A 2 e Jf m 
SS e^ep = e^ep and S$ ef 2 e^ = ef 2 e^\ for A, e AW A 2 e A^ (t) (A.27) 



The parameters ti2, ri 2 are defined in (IA.18I) and A/f(4.) are given in (12 .4ft . 

The notation A e A/^g) means that the excitation A can be equal to any value in A/^g) 

except the value 1 t (i) since it was chosen as the vacuum at first level. 

For an arbitrary N excitations number, we write 

$(n ab P) = S^(k a ,k b )$(P), (A.28) 

where the matrix S^ b \k aj k b ) has the same meaning as above but acts nontrivialy only on 
the V a ® V b vector space. It can be written using permutations and projectors: 

h) = - (^t + pAW + pM^ + pFMj + e _ ifea ^ t + + 



ab ' ab J ' ab ' ab ' ab ' ab 

+ (t ab I^ + r ab P^) + (Ulft* + r ab P^) (A.29) 
where we defined 

^ = EEW and ^' = EE^' £,£' = Ar t ,A^,A7 t ,A7,. (A.30) 

The matrix S^ b \k a , k b ) satisfies the Yang-Baxter equation (12.81) . 

Now we can write the periodic boundary conditions. Let CV = njvi-.-rijvjv-i be a circular 
permutation. The sites L + 1 and 1 being identified, we have the periodicity conditions: 

^qc n [ x ~~ e Q{N)L, ~A] = *§*[x,3] . (A.31) 

In terms of $(P) this yields the condition 

$(PC) = e ifc ^> L <l(P) . (A.32) 

If we choose P = C^ - "' with j = 1, . . ., N, we can derive a system of equations on the 
coefficients $(z<i) which is called the "auxiliary problem": 

t )j $(td) = e^ L $(td) with f) j = sW lJ ...S§sW...sV ltj , j = l,...,N (A.33) 
where we omitted the arguments of the S'-matrices, S^ 1 = 5^(fc a , 

The Yang-Baxter equation for the S matrix implies that [t)j , fjfc] = for all j, fc, so that the 
new Hamiltonians f)j can be simultaneously diagonalized: we do it in the folowing section. 
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A. 2 Auxiliary problem, level two 

In order to simplify the calculations, we make the change: S$ —> —S^ = £12- The 
eigenvalue problem to be solved reads 

S j+1J ...S Nj S lj ...S j _ hj <f) = Aj<f) . (A. 34) 

Again we use the coordinate Bethe ansatz. At this level we have n + m — 2 types of different 
excitations: e 2 ^, e n ^, e 2 ^, ...,e m ^. We choose as the reference state (pseudovacuum) : 

N 

# = n^ (A.35) 

fc=l 

The eigenvectors for this auxiliary problem are given by 

K 

^)[A] = Yl ^Min^^ (A.36) 

X6[1,JV] n=l 
p q n m— 1 

i=2 i=2 i=p+l i=q+l 




The ansatz for ^^[x, A] = ^q[x, A] can be written as, for Xq{\) < £q(2) < ■■■ < xq(k) 
and Q G &k, 

v Q [x,a]=j2®a(pq,p~ i ) n Map®) n f*i( x P(i)) n ^(^«) n 

P ieM^t ieMitt ieM^i ieM^ 

(A.37) 

where P = P^P^P^P^i is a factorized permutation such that P^, P^, -P74 and are 
the sets of permutations in &k permuting only particles of type tt j, tt T? ^ 4- an d 4- 
respectively, e.g. Pa(«) = i for z ^ .M/t and Pa(J) e -^/i for j £ with A = 7r t> ^ t> 
7r 4, 7f |. We recall that the sets Ai n ^, M.^ v ■Mn i and -^th are defined in section [31 

The eigenfunctions f x {\), g x {b), h x (a) and h x (a) correspond to one-particle solutions and 
are defined as 

/.(A) = if f~ A Sm , fcm + \ ) 5 ^) = & x ; ^(a) = ^(a) = 5(x-a). (A.38) 
^ V isinAm+i + zA- |y 

Although the functions /i and h are identical, they will be associated to different kinds 
of particles (tt I" and tt 4, respectively). Since they will lead to different Hamiltonian eigen- 
values (see below (1A.44P ). we have chosen to distinguish them in ( 1 A. 3 7ft . The form of the 
eigenfunctions (IA.38I) is supported by the fact that the eigenfunction in (1A.33P should be 
independent of the index j. 
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We can remark that the Bethe roots a\ for % G M.^ or M.^ corresponding to and 
particles are already quantized on the small chain. In order to have an independent number 
of eigenvectors flA.361) . one should take the following conditions into account: 



di < a i+ i for i G M., 



't(40 



(A.39) 



that leads to the fact that in (I A. 3711 . P^ and are fixed: Pn tU) Q(j) = const (see below). 

The action of the auxiliary Hamiltonian ( 1A.33I1 on the wavefunction can be calculated 
using the relations given in the previous section. In the recursive representation of the 
Hamiltonian (for details, see |18j ) 



A' 



• • • S^^iA] = H^ Q [x, A] J] e 

x6[l,JV] n=l 



A n Al) 



(A.40) 



one can derive the following recursive relations between the coefficients Hj ^q[k, A] for de- 
creasing k. Nontrivial relations occur when acts on ^q[x, A] for which some coordinate 
x p is equal to j and/or j — k. One gets two sets of relations. The first one reads: 



H^ Q [x',j;A] 



t^Hf-^qW, j;A}- r^H^^g^ ', j - k; A] for p G M.^ 



Hf-^Q^'.j-k^A] 



(fc-i) 



V 



v 



•Qi 



for p g M^i 
for p G Mrf 
for p G M ni 
(A.41) 



The above relations are invariant if we exchange j with j — k except the last one, that is 
transformed as follows: 



Xp 

I 



H ( ; ) y Q {x', J lk;A] = -e 
The second set of relations is 



•■i-kjj 



Xp 

I 



t, A] for p G 



(A.42) 



k;A] 



< Mk-l) 
j 



'^'[x'j -k,j;A], 



V Q >[x!,j - k,j;A}, for A pi ,A P2 G M nai M f 
^H^Q^JJ-k;!}-- - 

for pi G M^ui), P2 e Mnift) 
-e- ik i-*H<f " 1} ^q[x', j, j - fc; A], for Pl G p 2 G A^ m) 

-e^-Hf - 1) \P g [x', j, j - fc; A], for Pl G M^W- P2 e 



*q[x', j, j - fc; A], for Pl g jM^), p 2 e M H(t) 
-H^hqlx'JJ -k;A], for Pl G M ntii) , p 2 6 M^(t) 
" 1) *q[x / , j, j - A;; A], for Pl G -M^ (t ), P2 G -M na) 



(A.43) 
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In eqs. (1A.41I) - (1A.43|) . the vector x' corresponds to the vector x without the components x p 



or x pi ,x P2 and does not contain any position x q = j or j — k. 

Using these relations, we can then apply the whole product of S'-matrices flA.34j) on a one- 
excitation function. We get for x ^ j 

Hf^UW = o-i(A)/.(A) and flf ~ l) g x {b) = bg x (b) 

H<f- l) h x {a) = -h x (a) and Hf~ 1] h x {a) = -e~ lk *h x (a) (A.44) 

and for x = j 



N 



flf _1) /i(A)= II ^(A)/i(A) and Pf ~%(b) = b l ~ N g 3 (b) 

Hf-^hjia) = (-l^hjia) and Hf~%(a) = {-l^^e^^h^a) (A.45) 

where 



Now we will consider the iT-excitations eigenvector ( 1A.36I) with the ansatz ( 1A.37I) . There 
are three different possible cases: 



I) there exists a Bethe root a\ = j for some I G M.^, with j being the index in flA.331) . 
II) there exists another Bethe root a m = j for some m G A4 n ±, 
III) there is no such Bethe roots. 
We detail the calculations for each of these cases. 

I) There exists a Bethe root cii = j for some I G Ai^ 

In this first case, let us introduce a set of integers {c^ G [1, K]}i e x jrt such that we have 
Q(oti) G Ai n -[- The set is ordered: oti < aj+i. We also define oti such that XQ( ai ) = j. 
The Hamiltonian acting on the wave function gives the following result: 



X Q(";) 

x Q(<*i) 

_l\N-l-N Ue ik 3 N n 



rj-iV-lJ, T , / + s 

ti- W Q (x;a; (3 ( ai) ,...,iQ (ctI _ l) , j , x Q{ai+l) , ...,x Q{aNf j) 

4- 



'qc^Ix;xq( Q2 ), j ,XQ( Qi+l) ,XQ( Qi+2 ), ...,Xq (qi )J 

= A i*Q( x ; ^Q(ai), -(^(aw),;, ^Q(a i+1 )' -> X Q(*N t j) ( AA7 ) 

and finally we get the condition on $^(PQ, P _1 ): 



$ A (pg,p-!; 



A^-l)^- 1 -^^-^^ . (A.48) 



. . Ot jy. — 1 

In equality ( 1A.48I) . we have introduced CVt = 11^ ...Ua^ an d f° r P = ^rt^H^^-^i' 
P' = P^P^P n where P^ f = n^.JIJ^. 
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If we change the sector Q — > QH a b, which means that we interchange two identical 
particles xq^ and Xq^, the coefficient Q A (PQU a b, P -1 ) should remain invariant: 

$ A (P'QTl ab ,P'- 1 ) = $ A (PQ,P~ 1 ) with Q(a),Q(b)eM^ (A.49) 

where P' = P^P^P ^P^ and is such that P^ = P^U Q{a)Q(b) . 
The relation (j A. 48ft can be simplified to: 

with CVf = njy- ...n^' r \ Now, taking the product of this equation by itself N^ n times and 
changing P C^P -)■ C^ t P -)•...-»■ C^P, we obtain 

Aj (-1)^-^-%^* = e 2 "^ (m t7r = 1, AV) (A.51) 



It determines the eigenvalue of the auxiliary problem for spin up 7r-particles. We also recall 
that the result has to be multiplied by (— l)^ -1 in order to come back to the initial problem 
due to the sign change done at the beginning of the section, see (IA.34[) . 

II) There exists another Bethe root a m = j for some m G Ai n ± 
Following the same steps as in the first case, we obtain 

$ A (PQ,P-i) ~ A ' [ b ] 

where P = P^P n P n , P' = P^P^P n P n and C nl = n^il^ -1 . 
Finally using the same trick for the product of equations, we obtain 



A ^_^ N - Nue -ik AN ^ +Nim) = ^i N - ^ ^ = ^ ^ N ^ (A ^ 53) 

It allows to determine the eigenvalue of the auxiliary problem for n 4, particles. Again, 
this result should be multiplied by a factor (— l)^" 1 to come back to the initial problem. 

For identical particles we have the relation: 

^(P'QIU P'" 1 ) = $ A (PQ, P- 1 ), with Q(a),Q(b) e (A.54) 
where P' = P^P^P^P n and P^ is such that P^ = P n iYl Q ^ )Q{b) . 

Ill) There is no such Bethe root 

In this last case, we have 7r-particles with spin up or down distributed anywhere on the 
chain but the site j. Acting with the Hamiltonian ( 1A.33I) on the wavefunction ( 1A.37I) we get 



•Ik 



Hf- i) * Q ix,A}=(-i) N ^ +N ^ n e -*-« n aj(Xi) n &** q m], if x^a.ss) 

l&Mni i&M^ i&Mai 
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with the following conditions on the coefficients <& A (PQ, P 1 ): 
f) for all PQ{i),PQ{i + f ) in M n 

§ A {PQH ii+ i, (PIIq^q^+i))- 1 ) _ iXpQfj) - iXp Q (i + i) - | _ _ ± 

®a(PQ, P- 1 ) = iXpQii) ~ ^Ap Q(i+1) + § " a «J«)«J(i+i) 

2) for all PQ(i),PQ(i + 1) in we impose 



$ a {pqii u+1 , {pu q 



-1- 



^(pg,p-i) 

3) for all PQ(z) G and PQ(z + 1) G Al n 



(A.57) 



$^(PQ,P-i) (A ' 58j 

the same relation holds for all PQ{i) G Af^f, PQ{i + 1) G Al^ and changing Q — >■ Qn ii+ i. 

4) for all PQ(z) G and PQ(z + 1) G M n 

<S> A (PQT[ ii+1 ,p- 1 ) -ik a , \ /. , Q x 

$^(PQ,p-i) ~ PQWfr ap QW l A PQ(i+i)J ( A - 59 ) 

5) for all PQ(i) G Al^ and PQ(i + 1) G M n 

$T7pg^pTn = e PQW ^a PQ(l) (A PQ(m) ) (A.60) 

6) for all PQ(i) G Al^ and PQ(z + 1) G M n 



<$> A {PQ,P^) 
7) for all PQ(i) G Al^ and PQ(z + 1) G Al n 

^(pqil^p- 1 ) 
$^(pg,p-i) 



b PQ(i+i) ( A -6i; 



5pQ(i+i) (A.62) 



In addition, for all PQ(i),PQ(i + l) in M.^ or Al^ we impose the condition for identical 
particles: 



d^(PQ, (PU Q 



-1 - 



<S> A (PQ,P^) 

For the perodicity conditions, there are again two different subcases: 
III-A) a 7f ^-particle is on the j site, 
III-B) a 7f f-particle is on the j site. 
Each case leading to slightly different conditions, we treat them separately. 



(A.63) 
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III-A) A 7r J,-particle is on the j site. The periodic boundary condition on the coeffi- 
cients <&a{PQi P _1 ) comes from the relations 

tff-^Qfx, j;A] = (-l)^t+^ "Q e -ik ai "Q a . {Xi) "Q fc.^).^ andx/j 

ie-M^ ieM*t jeA^ 

Using relations (IA.560 to (1A.63|) we can calculate _1 ^\I/q[x, A] to get the condition 

*iTpQ K f^? = [bpQ(K)]N with cf *= n k- n £" 1 - 

We start with Q = such that 

~ /«, fori^jM n / m = P(i) e [1, IT] D A% i ^ 

[j e M n , for % e .M 5 ^P^Q^ 1 ) ee P{%) e M n , % e M n 

(A.65) 

and 

pnr _ n P(i) n P(2) u p( - K ~ N ^u p( - K ~ N ^ +1) n p{K - l) pn (A rr\ 

PQC K - Ll p{K) lL p{Ky .Al p{K) ll p{K) -ttp {K) PQ- (A.bb) 

Using relations (1A.61[) . (1A.62j) and "conjugated" (1A.58|) when Q — > QHa + i, the coefficient 
^(PQC^P- 1 ) simplifies to: 



= [^wl^^^^^n^J-^^^.n^-^pg, p- 1 ) . (A.67) 

Finally we have the equation 

*a(K J^ +1) ...nJ ^pq, p- 1 ) = \b Pw ] N -^ +N «+ N ^9 A (PQ, p-i) (A . 68) 

In the next step of the procedure, we introduce the vector ^^^(P): 



*MP) = E E M^tf) II e i (A.69) 
where G = span{(q + 1) J,, (m — 1) |}. The particle order is chosen such that for Q' = 



id, the vector ILeM^ ef* is equal to e (q+1)4 - <g> ... <g> e (c,+1)J - <g>... <g> e (m_1)i (g) ... <g> e (m ~ 1)J \ The 
permutation Q' is the image of Q^i when one shifts by K — the space of integers on 
which the permutation group & m _ acts. 

The relations flA.57j) and flA.63j) can be gathered in $^(P): 

$ n (n ab P) = Pa6^(P) (A.70) 

where V a b is the permutation acting on the particles located at the a and b positions in 
$7t[_(P). The periodicity problem can then be rewritten in the following form 

^(^"^•■■n^P) = [b P(K) } N - K+N ^n(P) (A.71) 
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which can be simplified to 

Pp {K -N i9+ l)P(K)...Pp(K-l)P(KAl(P) = [b P{K) ] N - K+N ^$ n (P) (A.72) 

Choosing P = C^™ with CV; = H% Nl " +1 ...U^ -1 , we get the following Bethe equation 



(b m ) N - K+N ^ = exp [ ^1 n,\ ,me M n (A.73) 

with 1 < m < ... < n N ^_ Num _ 1} < (A. 74) 

In order to get all quantum numbers which characterize the eigenfunction and obtain the 
right number of states, one should solve the permutation problem. We will take the result 
from [T3] (section on permutation problem of (g/(n|m) © <?/(2)) model). We get additional 
sets of integers 

1 < n?> < n { 2 k) < ... < nSi q+1)+ ...+^ fc _ 1 , < N liq+1) + ... + N ik , k = q + 2,...,m-2. (A.75) 

One can remark that if we set k = m — 1 in the above condition, we recover the relation 
0A.74p . This result ends the first subcase. 

III-B) A 7r ^-particle is on the j site. The periodic boundary conditions on the coeffi- 
cients $^(P(5,-P _1 ) can be written analogously to the previous case 

®$ P ( PQ K p P iP = H<n(*PQ(K)) fOT Q such that Q{K) e M n . (A.76) 

We proceed as in the case III-A, but instead we choose Q = Q^C^- u tir . Here 
Qt^ defined similarly to (5*4.: Qw[{i) G -M^t if i G and Q#t(i) = * f° r ^ ne remaining 
indices. The cyclic permutation is chosen such that all tt and 7f down particles are moved 
to the begining from their original ordering for Q = id, that is C K u u {i) = (i + N^ n + 
A^jf) mod x for any z. This choice implies in particular the following relation: 




pqc k = ( n o( n n s;u,,)( n n 



^(0 

P(AT t7r +Af t# )y ^ 11 ^P(N U +N^)J \ 11 "P(JV t7r +jV ti ) 
iGA4 n «G.M n i€M^ 

jP(N u +l) P(N tn +N^-l) 
■■ P(N U +N n ) 

p n K-N u -N^ 



X 



where P = P^Q^f and PQ = P^PC K u u . Here, due to the choice of Q, P^ and 
P^ are equal to identity. The symbol ^P{N r +JV t -) denotes the ordered product of 
permutation: ...n*-.„ r ar ,n*i», , T 

r P(N U +N n ) P(N tn +N n ) 

Next, we use relations f TA~58|) . (1AT59]) and f[Q0|) to simplify the coefficient ^ P p Q ? K : 



p- 



ZaIPQCk^- 1 ) = n e^a a XX P{Nu+Nn) ) f[ bj l J] e"^a ai (A P „ 

r P(JV tw +l) P(N u +N n -l) vn D _i 
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and finally we get 

TT J k a, „ ( \ _ \ TT TT ~ f i _ 



n e ik ^a a xxp (Nu+Nn) ) n ^ n & ^^(v* 



N 



Now we introduce similarly as in previous case the vector $^(P): 



**t(^)=£ e ^(^qo n e ^ Q ' w ( a - 8 °) 



where F = span{(p + 1) |, n t}- 

The relations flA.561) and flA.631) can be gathered in $^(P): 

$(n ab P) = a- b 1 V a MP), with P- I (a)-P" 1 (6) = -1 (A.81) 

where V a b is the permutation acting on the particles situated on a and b positions in $^(P). 
Therefore, the periodicity problem can be rewritten as follows: 

N 

X m V P{N ^ +1)m ...Vp {Nu+Nn _ x)m ^{P) = ;Q^(A m )<Mn (A.82) 



i=i 

with m = P(N U + Nfa) and 

X m = J] a ml J] e^a ai (A m ) J] 67 1 J] e-^c^AJ (A.83) 

Choosing P = C^ n with C^f = II^^^ ...1T^^^ 1 , we find the Bethe equations: 

J] a mi II e ^(A m ) J[ K 1 J] e"^ a ai ( A m ) e "t* ^ = JJ ^ (A,„) 

with m G A^fs- and 1 < fi\ < ... < n Ntii - Ntn < (A. 84) 

The problem of the complete characterization of the eigenfunction and of state counting 
arises in a similar way as in case III-A. We get additional sets of integers: 

1 < n[ k) < n { 2 k) < ... < n^ (p+i)+ ... +JVt(fc _ i) < iV t{p+1) + ... + iV> fc , k = p + 2,...,n-l. (A.85) 

Again, we recover relation f!A.84j) if we take k = n in the above condition. 
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